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In  vision  of  everyday  scenes^ features  requiring  detection  are  frequently 
observed  in  the  presence  of  suprathreshold  background  structures. 
Detection  of  such  features  is  a  contrast  discrimination  task  and  is  often 
necessary  for  the  subsequent  process  of  recognition.  In  order  to  provide 
a  description  of  this  task,  contrast  discrimination  measurements  were 
determined  for  targets  with  luminance  profiles  which  were  localised  in 
both  space  and  spatial  frequency.  The  investigation  extends  earlier 
work  on  this  topic  by  measurement  of  contrast  discrimination  levels  for 
different  base  contrasts,  sites,  luminances, and  aspect  ratios  of  the 
targets.  A  model  is  proposed  to  describe  the  contrast  discrimination 
process,  and  an  example  is  given  of  a  simple  application  of  the  model  to 
the  determination  of  the  number  of  discrimlnable  steps  in  contrast  as  a 
function  of  the  spatial  frequency  of  a  sinusoidal  grating  target.^ 

Approved  for  issue: 

G  H  Cockett 

Head  of  Applied  Physics  Department 

Reprint  with  minor  amendments  of  a  paper  published  in  Biological  Cyber¬ 
netics,  Vol  34,  1981. 
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1  INTRODUCTION 

This  report  is  the  first  in  a  planned  series  concerning  work  which  is 
primarily  directed  at  describing  visual  performance  with  imaging  systems. 


Much  research  in  spatial  vision  has  been  concerned  with  the  measurement 
and  modelling  of  performance  levels  for  the  detection  of  structures 
presented  on  uniform  background  fields.  On  the  basis  of  the  lead  given 
by  Schade  (Ref  1),  many  experiments  have  been  designed  around  a  linear 
systems  analogy,  and  as  such  have  used  sinusoidal  gratings  as  stimuli  in 
order  to  define  some  assumed  fundamental  responses  of  spatial  vision 
(Refs  2,3).  The  signal  transfer  function  of  the  system  is  taken  to  be 
linear  for  small  signal  amplitudes.  At  a  given  ambient  light  level,  it 
is  assumed  that  the  visual  system  can  be  represented  by  a  single  spatial 
filter:  a  single-channel  model.  However,  there  is  considerable  evidence 
that  several  mechanisms  may  mediate  threshold  detection  leading  to  the 
formulation  of  multiple-channel  models  (Refs  4-7).  In  addition,  the 
assumption  that  the  contrast  sensitivity  function,  as  derived  from  detec¬ 
tion  threshold  levels  of  sinusoidal  gratings,  represents  a  modulation 
transfer  function  in  the  conventional  sense  has  been  questioned  (Refs  8- 
10).  Such  studies  show  that  the  band-pass  attribute  of  the  spatial 
frequency  response  function  becomes  progressively  less  frequency  selective 
as  contrast  Increases  (Refs  11-13),  as  if  apparent  contrast  is  equal  to 
physical  contrast  minus  contrast  threshold.  A  process  of  neural  'deblurring* 
has  been  proposed  by  Georgeson  &  Sullivan  (Ref  8)  to  explain  this  effect. 
Snyder  and  Srlnlvasan  (Ref  14)  consider  a  more  general  Interpretation  in 
which  the  visual  system  compensates  for  both  optical  degradation  and 
noibe  contamination  in  the  signal  to  noise  discrimination  task.  A  non¬ 
linear  signal  transfer  function  of  the  decelerating  type  (eg  logarithmic) 
will  cause  the  apparent  contrasts  of  gratings  of  all  spatial  frequencies 
to  tend  towards  equality  (Ref  15),  though  this  behaviour  is  not  consistent 
with  the  contrast  linearity  observed  by  Cannon  (Ref  10). 

Suprathreshold  contrast  measurements  seem  to  be  at  variance  with  data 
from  single-cell  recordings  which  can  show  spatial  frequency  dependent 
(suprathreshold)  responses  at  most  stages  throughout  the  primary  visual 
pathway  (Refs  16-18).  This  may  mean  that  significant  spatial  processing 
occurs  at  levels  more  central  than  those  as  yet  examined  by  single-cell 
recording,  that  populations  of  cells  presumably  active  in  psycho-physical 
experiments  have  mass  properties  which  are  not  readily  derived  from  single¬ 
cell  responses,  or  possibly  that  some  cell  response  other  than  spike 
frequency  is  the  critical  parameter  (eg  the  statistics  of  the  relationship 
between  individual  spikes  (Ref  19). 

In  spite  of  the  evidence  which  reveals  flaws  in  the  general  applicability 
of  single-channel  models  of  spatial  vision,  such  models  have  been  widely 
used  in  the  design  of  Imaging  devices,  especially  electro-optic  systems. 

This  is  presumably  because  such  models  capture  the  essence  of  spatial 
vision  in  that  they  provide  a  decreasing  response  as  spatial  frequency 
Increases  (which  implies  an  acuity  limit),  and  give  a  poor  response  to 
diffuse  targets.  In  addition  the  models  have  a  simple  construction  and 
are  fast  in  application. 
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Models  consisting  of  several  channels  are  necessary  in  order  to  embrace 
the  effects  on  detection  threshold  levels  of  subthreshold  structures 
(Refs  6,7)  and  masking  background  fields  (Refs  20,  21).  Such  models  are 
consistent  with  much  experimental  evidence  (Refs  22-25).  In  vision  of 
everyday  scenes,  however,  features  requiring  detection  are  often  observed 
in  the  presence  of  suprathreshold  background  structures.  The  detection 
of  such  features  is  necessary  for  the  process  of  recognition.  Therefore 
a  description  of  suprathreshold  vision  is  not  only  useful  for  specifying 
the  appearance  of  objects  (cosmetic  quality)  but  also  crucial  for 
determining  those  features  of  an  object  which  are  available  for  subsequent 
recognition  by  the  trained  higher  brain:  the  primal  sketches  of  Marr  (Ref 
26). 

Assuming  that  magnitude  estimation  (Ref  10)  and  contrast  halving  (Ref  9) 
provide  correct  descriptions  of  signal  transfer  functions  in  spatial 
vision,  it  seems  that  perceived  contrast  is  equal  to  objective  contrast 
minus  contrast  threshold.  If  in  addition  the  statistical  spread  of  supra¬ 
threshold  responses  is  known  for  a  fixed  target  contrast,  then  a  first 
step  can  be  made  towards  modelling  suprathreshold  spatial  vision.  Such  a 
model  can  be  applied  to  several  problems:  the  prediction  of  just-noticeable- 
differences  in  measures  of  image  quality  (eg  system  MTF) ;  the  accuracy  of 
visual  accommodation  required  to  avoid  loss  of  image  quality  (Ref  27); 
and  the  design  of  digital  image  coding,  transmission  and  displays  systems 
in  which  luminance  levels  are  quantised  and  interact  with  spatial  sampling 
intervals  (Refs  28-32). 

The  aim  of  the  present  investigation  is  to  measure  contrast  discrimination 
levels  for  different  contrasts,  etc,  sizes,  luminances  and  aspect  ratios 
of  targets  with  triphasic  spatial  profiles  (Fig  la).  Such  targets  are 
localised  in  both  space  and  spatial  frequency,  so  avoiding  probability 
summation  effects  (Refs  25,  33,  34)  while  retaining  moderate  spatial 
frequency  selectivity.  The  investigation  extends  the  results  of  Nachmlas 
&  Sansbury  (Ref  35),  Kohayakawa  (Ref  36)  and  Kullkowskl  (Ref  9).  A 
further  aim  was  to  place  on  a  firm  experimental  basis  assumptions  made 
by  Mezrich,  Carlson  &  Cohen  (Ref  30)  in  their  calculations  of  the  theore¬ 
tical  number  of  grey  scale  steps  as  a  function  of  spatial  frequency,  and 
the  minimum  number  of  bits  required  for  faithful  image  reproduction  on 
displays. 

2  METHODS 

Targets  were  generated  by  a  Microconsultants  INTELLECT  framestore  and 
image  processing  system  controlled  by  a  Computer  Automation  LSI-2  computer. 
The  frame  store  has  a  resolution  of  512x512  picture  elements  (pixels) 
with  an  8-bit  grey  scale.  Output  from  the  framestore  drives  a  digital  to 
analogue  converter  at  50Mblt/sec  enabling  targets  to  be  displayed  on  a 
broadcast  quality  monitor  (FYE  TVT  Model  LDM15  fitted  with  a  green  (OH) 
phosphor;  chromaticity  coordinates  x-0.270,  y*0.602).  The  observer  viewed 
the  display  from  a  range  of  3.5m  with  his  head  resting  in  a  close  fitting 
support.  At  this  range  the  screen  subtended  4.7  deg  horizontal  and  3.2 
deg  vertical.  The  pixel  subtense  was  0.55mln  horizontal  and  0.375min 
vertical.  A  rectangular,  uniform  surround  field  (14.7deg  horizontal, 

10.8deg  vertical)  was  provided  by  means  of  mat,  white  card  and  a  projector 
with  an  interposed  filter  solution  made  from  cupric  chloride  and  potassium 
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produce  an  approximate  match  to  the  colour  and  luminance  of  the  monitor 
display.  The  monitor  luminance  was  set  at  28cd.m’'^.  In  order  to  maintain 
this  value  of  luminance  the  monitor  AC  supply  was  provided  by  a  Cintec' 
voltage  and  frequency  stabiliser.  lower  luminance  values  were  achieved 
by  placing  large  neutral  filters  (titanium  on  glass)  close  to  the  eyes  of 
the  observer.  The  experiments  required  the  presentation  of  contrasts 
ranging  from  about  0.001  to  0.8  (equation  1).  By  selecting  different 
settings  of  the  monitor  brightness  and  contrast  controls,  while  maintaining 
constant  mean  luminance,  three  overlapping  contrast  scales  were  obtained. 
Thus  sufficient  luminance  resolution  was  retained  in  the  mapping  from  the 
frame  store  grey  scale.  Each  contrast  scale  was  calibrated  using  a  Gamma 
Scientific  telephotometer,  model  2009.  Appropriate  non-linear  corrections 
were  applied  so  that  the  contrast  as  specified  in  the  controlling  computer 
program  was  the  same  as  that  which  appeared  on  the  monitor. 

Most  values  of  contrast  discrimination  were  obtained  by  using  a  two- 
alternative  forced-choice  paradigm  designed  to  achieve  a  71%  probability 
of  correct  discrimination  between  two  targets  of  different  contrast. 

Coded  tones  indicated  correct  and  incorrect  responses.  (Methods  employing 
'same'  or  'different'  responses  were  tried  but  were  found  to  be  unreliable 
because  nominally  identical  targets,  other  than  point  sources,  almost 
always  appear  different  during  a  given  presentation  owing  to  random  point- 
to-point  fluctuations  in  visual  responses.)  Values  of  contrast  were 
presented  by  using  a  randomly  interleaved,  double  staircase  procedure 
(Ref  38)  with  adaptive  steps  (Ref  39)  so  that  the  region  of  71%  correct 
discrimination  was  achieved  rapidly  with  large  contrast  steps,  leading  to 
6mall  contrast  steps  as  the  threshold  value  became  more  localised.  A 
number  of  preliminary  presentations  (typically  40)  were  made,  until  the 
step  size  in  contrast  difference  for  both  staircases  fell  below  20%. 

After  this  point  40  further  presentations  occurred.  The  final  value  of 
contrast  discrimination  was  determined  by  averaging  the  40  values  of 
contrast  difference  for  both  staircases.  In  order  to  determine  the 
experimental  errors  appropriate  to  this  procedure,  some  double  staircase 
procedures  were  repeated  10  times.  Standard  deviations  of  contrast 
discrimination  were  obtained  by  this  method. 

The  main  experiments  used  side-by-side  presentation  of  the  two  targets. 

The  observer  was  instructed  to  direct  his  gaze  from  one  target  to  the 
other  during  the  2  second  presentation  and  ignore  on  and  off  transients. 

The  centre-to-centre  target  separation  was  0.92deg  for  all  triphasic 
targets  except  for  the  largest  target  used  (centre  bar  width  35min)  for 
which  the  separation  was  increased  to  1.84deg.  The  Interval  between 
target  presentations  was  at  least  3n  during  which  time  the  observer  viewed 
a  uniform  field.  A  small  dark  reference  crons  or  spot  was  provided  in 
the  centre  of  the  monitor  field  in  order  to  aid  target  location  for  near 
or  below  threshold  targets.  Slde-by-alde  presentation  of  targets  was 
used  in  preference  to  temporal  forced-choice  methods  since  the  former 
relies  less  on  memory  and  more  closely  resembles  the  type  of  comparison 
of  image  quality  often  employed  in  system  design.  However,  a  short 
experiment  was  performed  using  temporal  forced-choice  presentation:  first 
target  on  for  2s,  uniform  field  for  Is,  second  target  on  for  2s,  then  a 
delay  of  at  least  3s  before  presentation  of  the  next  pair  of  targets. 
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Most  experiments  used  trlphsslc  targets  consisting  of  a  central  bar 
flanked  by  two  bars  of  opposite  sign  (contrast  profiles  are  shorn  In  Fig 
la).  The  examples  In  Figure  la  show  positive  triphasic  profiles:  a 
central  bright  bar  flanked  by  dark  bars.  The  alternative  polarity  of 
target  (negative  triphasic)  was  also  used.  The  targets  are  more  localised 
than  the  grating  patches  employed  by  Graham,  Robson  and  Nachmias  (Ref  34) 
and  therefore  some  spatial  frequency  selectivity  Is  lost,  though  the 
targets  are  more  typical  of  the  structures  encountered  in  real  objects. 

The  continuous  lines  In  Figure  la  are  records  of  luminance  obtained 
directly  from  the  monitor  screen  using  a  Gamma  Scientific  photometric 
microscope,  model  700-10  with  75u  silt.  The  fine  structure  displayed  In 
the  scans  is  due  to  phospher  granularity.  The  open  symbols  describe  the 
contrast  values  calculated  In  the  experimental  computer  program  with  a 
small  scaling  correction  applicable  at  high  spatial  frequencies  (  30c/deg) 
to  account  for  the  roll-off  in  whole  system  response.  This  reduced 
response  is  evident  for  the  narrow  targets  shown  in  Figure  la.  Target 
width  for  triphasic  profiles  is  defined  as  the  width  of  the  central  bar, 
le  the  distance  between  the  two  zero  contrast  crossings.  Target  contrast 
Is  specified  by: 

C  -  Il-IqI/Io  (1) 

where  L  is  the  luminance  of  the  middle  of  the  central  bar  and  Iq  is  the 
background  luminance.  For  a  positive  triphasic  profile  of  unit  width  and 
contrast,  the  equation  describing  the  profile  as  a  function  of  the  spatial 
coordinate  0<|x| <1.305  Is: 

C(x)  -  l-l.829x-l. 073x2+l. 447x3  (2) 

The  modulus  of  the  Fourier  transform  of  the  function  for  a  target  of  size 
D“4.4min  and  contrast  1.0  is  shown  by  the  open  symbols  In  Figure  lb.  A 
similar  transformation  of  data  from  a  photometric  microscope  scan  of  the 
monitor  is  described  by  the  continuous  line  in  Figure  lb.  The  good  agree¬ 
ment  between  theoretical  and  experimental  data  in  Figures  la  and  lb 
provides  a  useful  check  on  equipment  reliability.  For  a  target  of  width 
Dmln  the  maximum  ^m  in  the  corresponding  spatial  frequency  spectrum  is: 

Fra  ■  0.50/D  c/mln  (3) 

Two  observers  participated  in  the  experiments:  BM  aged  43  and  GJB  aged 
30.  BM  viewed  the  display  binocularly  and  was  not  informed  of  the  aims 
and  procedures  of  the  experiments.  GJB,  the  author,  viewed  with  the 
right  aye  only,  the  left  eye  being  occluded.  No  artificial  pupils  or 
spectacle  corrections  were  employed. 

3  RESULTS 

All  measurements  of  contrast  discrimination  relate  to  the  observation  of 
one  target  of  contrast,  OAC,  which  can  be  chosen  with  some  probability  P 
as  being  of  higher  contrast  than  a  similar  target  of  base  contrast  C. 

Except  for  the  data  shown  in  Figure  9,  the  probability  of  a  correct 
contrast  choice  in  all  experiments  is  approximately  0.71.  The  conventional 
contrast  threshold  for  a  given  target  therefore  corresponds  to  the 
condition  00.  The  targets  used  were  always  positive,  triphasic  profiles 
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FIG  1(a)  Continuous  lines  show  photometric  microscope  scans  of  some 
triphasic  profiles  employed  In  the  experiments.  The  open  symbols  represent 
the  corresponding  values  as  defined  In  the  controlling  computer  program. 
Rach  target  has  a  nominal  contrast  of  1.0  and  the  number  adjacent  to  each 
curve  defines  the  particular  target  width,  min  of  arc,  as  specified  by 
the  width  of  the  central  (bright)  bar.  The  microscope  can  scan  an  angle 
of  only  13  min  of  arc:  therefore  the  flanks  of  the  two  larger  profiles 
are  not  plotted.  The  small  reduction  In  contrast  for  the  narrowest  profile 
Is  due  Co  a  roll-off  In  whole  system  MTF  at  high  spatial  frequencies. 

(b)  Modulus  of  the  Fourier  transforms  for  experimental  and  theore¬ 

tical  triphasic  profiles  of  width  4. Amin.  The  total  width  between  half 
maximum  points  is  1.6  octaves. 
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( Fig  la)  except  where  identified  as  negative,  triphaelc  profiles  or  as 
linear  sinusoidal  gratings.  All  targets  were  orientated  vertically  and 
extended  the  full  height  of  the  display  3,2deg,  except  where  the  effect 
of  different  aspect  ratios  was  examined.  The  mean  display  luminance  was 
28cd.m~2  for  most  experiments.  In  some  figures  error  bars  are  shown  which 
represent  Usd  (N>10). 


Conventional  contrast  threshold  values  for  both  observers  are  shown  in 
Figure  2,  values  are  determined  for  different  target  sices.  The  'band- 
pass'  type  of  behaviour  as  observed  for  similar  measurements  with  sinus¬ 
oidal  gratings  is  apparent,  as  is  a  'cut-off  in  the  region  of  O.Oldag 
(approximately  equivalent  to  a  grating  of  spatial  frequency  of  SOc/deg). 

These  results  lend  support  to  the  statement  made  earlier:  the  trlphaalc 
targets  have  a  sufficiently  selective  spatal  frequency  content  to  produce 
visual  responses  similar  to  those  of  sinusoidal  gratings  without  the 
associated  difficulties  produced  by  the  presence  of  multiple  cycles  (Ref  33). 
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FIG  2  Conventional  contrast  threshold  values  for  positive  triphasic 
profiles  of  different  widths.  Filled  and  open  symbols  show  data  for 
Observers  BM  and  GJB  respectively.  Luminance  28cd.m~^. 


Measurements  of  contrast  difference  AC  determined  for  different  values  of 
hnse  contrast  C  are  shown  In  Figures  3a,b&c.  A  comparison  of  different 
experimental  techniques  is  shown  in  Figure  )a,  different  polarities  of 
triphasic  profile  In  Figure  3b  and  different  target  types  (triphasic 
profile  and  grating)  in  Figure  3c.  In  all  cases,  the  contrast  difference 
AC  Initially  decreases  for  increasing  base  contrast  C,  and  finally 
Increases  In  a  way  which  exhibits  an  approximate  linearity  of  AC  against 
C:  the  contrast  Weber-Fechner  region.  Overall,  the  type  of  variation  of 
AC  with  C  Is  the  same  as  that  observed  by  Campbell  &  Kullkowski  (Ref  AO) 
and  by  Nacbmias  &  Sansbury  (Ref  35).  The  initial  decrease  in  contrast 
difference  AC  has  been  termed  subthreshold  addition  (Ref  6).  The  minimum 
value  of  AC  (Fig  3)  occurs  at  a  value  of  base  contrast  C  which  Is  in  the 
region  of  the  conventional  threshold  level  for  the  target. 
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PIC.  U  Contrast  dlscrlalnatlon  valuas  obtained  for  different  widths  of 
triphasic  targets,  (a)  Observer  BM,  (b)  Observer  C.TB.  Values  were 
obtained  by  the  method  of  elda-by-alde  forced-choice  using  positive  tri¬ 
phasic  profiles.  The  dashed  lines  ere  dre’m  with  a  slope  of  1.0  end 
contrast  Webar-Fachnar  fraction  of  0.08.  The  letter  value  is  an  average 
estimate  of  the  fraction  as  shown  In  Figure  7.  luminance  28cd.m-2. 
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The  target  of  width  0.15deg  was  sleeted  and  used  to  determine  the  depen¬ 
dence  of  contrast  discrimination  on  the  display  luminance.  Values  for 
one  observer  are  shown  in  Figure  S.  For  each  luminance  the  variation  of 
contrast  difference  AC  with  base  contrast  C  has  similar  characteristics 
to  those  already  observed  for  the  influence  of  target  width  (Fig  4). 

Thus  each  set  of  symbols  mvoes  along  a  45deg  line  so  that  the  placement 
of  the  set  depends  only  on  the  conventional  threshold  level. 


FIG  5  The  Influence  of  luminance  on  contrast  discrimination  of  positive 
triphasic  targets  using  the  side-by-side  forced-choice  method.  Velues 
were  obtained  for  luminances  of  28,  0.31,  C.032  and  0.0036cd.m“2  and  are 
denoted  by  open  squares,  filled  squares,  asterisks  and  crosses  respectively. 
Observer  GJB,  target  width  O.lSdeg.  The  broken  line  is  drawn  with  a 
slope  of  1.0  and  a  contrast  Weber-Fechner  fraction  of  0.08. 


It  la  clear  from  the  results  shown  in  Figures  4&5  that,  to  a  first  approxi¬ 
mation,  values  of  AC  in  the  contrast  Weber-Fechner  region  are  Independent 
of  the  conventional  threshold  values  when  those  values  are  altered  by 
changing  either  the  target  width  or  the  mean  luminance. 

The  effect  A  different  aspect  ratios  of  the  target  (length/width)  was 
briefly  examined  by  determining  the  minimum  detectable  contrast  difference 
for  a  target  of  width  0.15deg  and  contrast  0.8  as  aspect  ratio  varied 
from  the  maximum  attainable  for  this  target  width,  21.3  down  to  1.0. 

Values  for  both  observers  are  shown  in  Figure  6.  There  is  no  systematic 
variation  of  contrast  difference  AC  with  aspect  ratio. 
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FIG  <>  Illustrates  chat  Che  length  of  a  triphasic  profile  does  not  signi¬ 
ficantly  affect  the  level  of  contrast  discriminate  in  when  the  base  contrast 
is  within  the  Weber-Fechner  region.  Filled  symbols  and  open  symbols  show 
data  for  Observers  BM  and  GJB  respectively.  Base  contrast  0.8,  target 
width  0.15deg  and  luminance  28cd.m“2.  Method  of  side-by-side  forced-choice 
using  positive  triphasic  targets. 


It  has  been  t.hown  that  the  contrast  different  AC  is  approximately  propor¬ 
tional.  to  the  base  contrast  C  for  values  somewhat  greater  than  the  conven¬ 
tional.  threshold  level  (Figs  2-5).  The  constant  of  proportionality  W  in 
the  relationship  AOW.C.  Is  defined  as  the  contrast  Weber-Pechner  fraction. 
Values  of  W  were  estimated  from  the  data  shown  in  Figures  4a4b  and  are 
plotted  against  target  width  In  Figure  7  for  both  observers  (square 
symbols).  Also  shown  Is  the  value  for  negative  triphasic  targets  as  esti¬ 
mated  from  the  data  shown  in  Figure  3b.  It  Is  apparent  that  the  contrast 
Veber-Feclmer  fraction  W  varies  slowly  with  target  size.  Nevertheless,  a 
mean  value  of  W,  0.08,  was  calculated  by  averaging  over  the  data  for  both 
observers  and  all  widths,  because  such  a  ' sing 1 --number'  provides  a  useful 
designer's  rule-of-thumb  and  Is  a  sufficiently  accurate  representation 
for  many  practical  situations.  A  value  of  W  was  also  estimated  from  the 
data  given  by  Nachmlae  &  Sansbury  (Ref  35).  This  single  value,  determined 
for  a  grating  of  spatial  frequency  3c/deg,  Is  plotted  In  terms  of  the 
equivalent  trlpheslc  bar  width.  Since  the  experimental  technique  employed 
by  Nachmlae  &  Sansbury  Is  similar  to  that  used  in  the  present  Investigation, 
the  magnitude  of  the  difference  between  this  single  value  and  the  values 
obtained  In  this  Investigation  la  sufficiently  large  to  justify  the 
averaging  process  used  to  derive  a  mean,  practical  value  for  the  Weber- 
Fechner  fraction.  It  Is  also  of  practical  significance  to  derive  an 
average  curve  for  the  relationship  between  the  contrast  difference  AC  and 
the  base  contrast  C.  This  was  achieved  by  taking  data  for  a  target  size 
of  0.15deg  at  a  luminance  of  28cd.m~2  as  the  standard.  Data  for  other 
conditions  were  then  adjusted  by  eye  to  produce  a  best  fit  to  the  standard. 
This  adjustment  was  made  by  moving  each  data  set  along  a  45deg  line  so 
that  both  contrast  and  contrast  difference  are  scaled  by  the  same  factor. 

The  results  of  this  process  are  shown  for  the  two  observers  in  Figures 
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8a&b.  These  scatter  plots  were  then  averaged  In  bins  of  width  0.2  lu 
(log  unit)  along  the  abscissa  to  yield  the  square  symbols  shown  in  Figure 
8c.  The  crosses  and  dashed  line  in  Figure  8c  represent  a  theoretical 
function  as  described  in  the  Discussion  (Section  4). 


■  oil  bmB  pos,t,ve  polarity 

*  OBS  GJB  NEGATIVE  POLARITY 


7  Estimates  of  the  Weber—Pechner  fraction  obtained  from  contrast 
discrimination  data  for  base  contrasts  in  excess  of  the  conventional 
threshold  value.  Filled  and  open  squares  are  values  for  Observers  BM 
(data  from  Fig  4a)  and  GJB  (data  for  Fig  4b)  respectively.  The  fraction 
obtained  for  negative  triphasic  targets  (Fig  3b)  is  indicated  by  the 
asterisk.  The  cross  shows  the  estimate  obtained  from  data  given  by 
Nachmlas  and  Sansbury  (Ref  35):  it  is  plotted  against  the  equivalent 
triphasic  bar  width.  A  mean  value  for  all  target  widths  and  observers  is 
indicated  by  the  dashed  line  at  a  value  of  the  contrast  Weber-Pechner 
fraction  equal  to  0.08. 


Up  to  this  point  the  values  of  contrast  difference  have  been  determined 
at  only  the  71%  level  of  discrimination.  For  purposes  of  theoretical 
analysis  and  practical  application  it  is  crucial  to  determine  how  the 
probability  of  discrimination  varies  with  values  of  contrast  difference. 
For  conventional  threshold  values  the  psychometric  function  is  approxi¬ 
mately  constant  when  plotted  on  a  normalised  scale  in  which  the  contrasts 
are  divided  by  the  contrast  oppropriate  to  some  fixed  probability  (Refs 
41,  42).  There  are  only  minor  differences  between  this  representation 
and  the  one  in  which  the  scale  is  represented  by  log  contrast  (Ref  43). 
For  the  contrast  discrimination  experimenta  it  is  not  immediately  obvious 
whether  the  scale  variable  should  be  the  contrast,  OfAC,  of  the  adjust¬ 
able  stimulus,  the  contrast  difference  AC,  or  some  other  function  of  C 
and  AC.  In  order  to  identify  the  appropriate  transformation,  psycho¬ 
metric  functions  were  determined  for  a  triphasic  target  of  fixed  width 
0.13deg  and  luminance  28cd.m~2,  Three  values  of  base  contrast  C  were 
chosen:  OO,  0.03  and  0.2.  The  first  value  corresponds  to  the  conven¬ 
tional  threshold,  the  second  is  just  within  the  contrast  Weber-Fechner 
region  and  the  third  is  towards  the  end  of  this  region  (Fig  3). 
Probability  values  are  shown  in  Figure  9.  They  were  estimated  from  an 
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FIG  8  Combined  data  of  contrast  discrimination  for  triphasic  targets 
versus  base  contrast  obtained  by  scaling  both  coordinates  by  the  same 
amount  in  order  to  produce  a  satisfactory  fit  to  data  for  target  width 
O.lSdeg  and  luminance  28cd.m*2. 

(a)  Observer  BM. 

0>)  Observer  GJB.  The  combined  values  are  Indicated  in  the  legend 
and  were  obtained  from  Figs  4a&b  &  S. 

(c)  Squares  show  averaged  contrast  discrimination  data  obtained  from 
(a)  and  (b).  Filled  and  open  squares  represent  mean  values  for  Observers 
BM  and  GJB  respectively.  The  crosees  and  broken  line  give  the  theoretical 
function  as  described  in  the  Annex.  This  function  was  fitted  to  the 
experimental  data  by  setting  tha  conventional  contrast  threshold  level 
equal  to  0.006,  K}"0.1  and  K2«8.  The  value  of  Kj  produces  a  contrast 
Weber-Fechner  fraction  of  0.08.  The  asterisks  show  sequential  calcula¬ 
tions  obtained  as  K}  was  held  constant  at  0.1  and  K2  was  increased 
from  aero  until  satisfactory  agreement  occurred  with  the  experimental 
data. 
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FIG  9  Psychometric  functions  obtained  by  determining  the  probability  of 
contrast  discrimination  aa  a  function  of  contrast  difference.  The  method 
of  alde-by-elde  forced-choice  was  used  for  positive,  triphasic  targets. 
Values  are  shown  for  three  base  contrasts:  00,  0.03  and  0.2.  Filled 
circles  and  open  circles  show  data  for  Observers  BM  and  GJB  respectively. 
Each  value  of  probability  Is  based  on  approximately  40  presentations  per 
point.  Continuous  lines  show  theoretical  calculations  obtained  by  the 
method  described  in  the  Annex.  The  arrow  on  each  abscissa  Indicates  the 
contrast  difference  required  for  a  probability  of  contrast  discrimination 
equal  to  7 IX.  Target  width  0.15deg,  luminance  28cd.m“2. 


average  of  about  40  presentations  per  point.  The  continuous  lines  are 
calculated  from  the  theoretical  derivation  given  In  the  Annex.  The  arrows 
on  the  abscissa  denote  values  of  contrast  difference  giving  probabilities 
of  0.71.  Note  that  the  psychometric  function  for  aero  base  contrast  is 
steeper  than  those  for  base  contrasts  within  the  contrast  Weber-Fechner 
region. 

4  DI8CU88ION 

Contrast  discrimination  has  been  measured  for  triphasic  targets  of 
different  widths  (Fig  4),  polarities  (Fig  3b) ,  aspect  ratios  (Fig  6)  and 
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luminances  (Fig  5).  Overall,  the  results  can  be  efficiently  described  by 
a  single  curve  of  contrast  difference  AC  against  base  contrast  C  with  the 
condition  that  both  AC  and  C  are  scaled  by  a  constant,  determined  only  by 
the  conventional  threshold  levels  The  valve  of  conventional  threshold  In 
particular  can  be  manipulated  either  by  changing  the  target  width  or  by 
changing  the  luminance.  The  standard  curve  (Fig  8c)  can  be  fitted  approxi¬ 
mately  to  all  sets  of  data  of  AC  against  C  by  the  choice  of  an  appropriate 
scaling  constant  which  Is  applied  to  both  axes.  This  description  of 
contrast  discrimination  is  sufficient  .y  accurate  for  many  practical 
purposes,  but  it  must  be  remembered  that  it  is  not  precise  since  repetitive 
experimentation  can  yield  variations  in  the  contrast  Weber-Fechner  fraction 
with  different  target  widths  (Fig  7). 

For  base  contrasts  C  somewhat  less  than  the  conventional  threshold  level 
Go  (Fig  8c),  the  results  are  accurately  described  by  the  equation: 

C+AC  -  C^;  C/Co<0.6  (A) 

The  linear  relationship  between  AC  and  C  in  the  region  of  subthreshold 
addition  confirms  previous  Investigations  (RefB  6,  4A)  and  follows 
logically  if  the  existence  of  both  small-signal  linearity  and  a  neural 
threshold  device  are  assumed.  If  no  such  device  exists,  then  a  non¬ 
linearity  in  the  form  of  an  accelerating  function  (eg  a  quadratic  depen¬ 
dence  of  output  on  input  contrast)  is  postulated  (Ref  35)  in  order  to 
provide  values  of  AC  which  decrease  with  increasing  (small)  values  of  C. 

For  the  purposes  of  the  present  analysis  the  existence  of  neural  threshold 
devices  is  assumed  since  such  mechanisms  provide  a  parsimonious  descrip¬ 
tion  of  the  experimental  results  of  subthreshold  addition. 

For  base  contrasts  C  somewhat  greater  than  the  conventional  threshold 
level  (o,  contrast  discrimination  follows  a  Weber-Fechner  relationship  in 
the  contrast  domain: 

AC-W.C;  C/Cq>3  (5) 

An  equation  of  this  form  can  arise,  for  example,  from  a  logarithmic  trans¬ 
ducer  function  with  a  constant  noise  level  added  after  the  non-linearity. 
However,  perceived  contrast  seems  to  be  equal  to  actual  contrast  minus 
contrast  threshold  (Refs  9,  10).  Even  though  this  (linear)  transducer 
fucntlon  could  be  the  result  of  cascading  two  or  more  non-linear  functions 
(Ref  45),  the  simpler  assumption  of  linearity  at  each  atage  has  been  made 
in  the  present  analysis.  In  order  to  reproduce  the  contrast  Weber-Fechner 
relationship  of  equation  5,  it  is  therefor”  necessary  to  Include  a  noise 
process  with  a  magnitude  which  is  proportional  to  the  contrast  (signal) 
level.  In  addition  another  it  Use  process  must  be  Included  which  has  a 
level  proportional  to  the  ccaventional  threshold  level  (<>;  ie  a  constant 
noise  level  for  a  given  target.  This  further  noise  source  ie  needed  in 
order  to  prevent  the  derived  psychometric  functions  from  becoming 
extremely  steep  at  low  base  contrasts,  especially  when  plotted  against 
log  contrast  difference  log  AC.  The  two  noise  sources  are  assumed  to  be 
uncorrelated  and  have  Gaussian  probability  distributions  (Fig  10b).  The 
noise  levels  can  therefore  be  described  by  equivalent  standard  deviations 
o.  For  the  first  noise  process  set  the  noise  level  Oj  by: 
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FIG  10  The  suprathreshold  contrast  model  used  to  derive  theoretical 
estimates  of  contrast  discrimination  (see  Annex). 

(a)  Shows  the  schematic  modelling  Illustrating  Interactions  of  sigral, 
noise  and  the  threshold  device. 

(b)  Probabilistic  responses  to  target  contrasts  C  and  C+AC.  The  dashed 
line  indicates  the  contrast  level  set  by  the  threshold  device.  All 
response  values  are  related  to  equivalent  input  contrast  levels. 

and  for  the  second  noise  process  set: 

°2*VK2*-cO  <7> 

where  Ki  and  K2  are  approximately  independent  of  stimulus  parameters  such 
as  slae  and  luminance.  The  conventional  threshold  contrast  is  Cq.  The 
combined  noise  level  from  the  two  uncorrelated  sources  is  then: 

-  I^.OCj.Cq^+C2)*, 

and  when  scaled  with  respect  to  Cq  is: 
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ST-  Kr(K2+T2)^  (8) 

where  S<j>  -  a/  Cq 

and  T  -  C/Cq 

A  model  was  devised  both  to  describe  the  contrast  discrimination  data  and 
to  be  consistent  with  subthreshold  summation  data.  It  is  shown  schematically 
in  Figure  10a.  Many  such  arrangements  can  be  used  within  the  components 
of  multiple-channel  models  (Refs  4,  25)  and  for  separate  processing  of 
different  contrast  polarities  (Refs  46,  47).  Such  parallel  processing  is 
not  described  specifically  in  the  present  analysis:  it  is  assumed  only 
that  some  unspecified  process  defines  the  conventional  threshold  level 
Cq.  An  equation  is  derived  in  the  Annex  which  gives  the  probability  of 
discriminating  targets  of  contrasts  C  and  04-AC  in  terms  of  the  conventional 
threshold  level  Cq  and  the  constants  Kj  and  K2«  The  average  val  le  of  the 
contrast  Weber-Fechner  fraction  W  in  equation  5  is  taken  as  0.08  (Fig  7). 

This  fixes  the  value  of  at  0.1  (Annex).  The  value  of  K2  was  then 
adjusted  stepwise  until  the  theoretical  variation  of  AC  with  C  gave  satis¬ 
factory  agreement  with  the  experimental  data  (Fig  8c).  The  value  chosen 
for  K2  is  8.  Note  that  If  only  one  parameter,  Kj,  is  used  (ie  K2“0) ,  the 
maximum  enhancement  in  sensitivity  over  the  conventional  threshold  level 
(Fig  8c)  is  over  1  lu,  whereas  the  experimental  data  indicates  about  0.6 
lu.  In  addition  and  as  noted  earlier,  very  steep  psychometric  functions 
occur  as  the  base  contrast  C  becomes  small.  The  complete  theoretical 
variation  of  AC  with  C  is  shown  by  the  crosses  and  dashed  line  in  Figure 
8c.  As  the  base  contrast  becomes  small  the  discrimination  process  tends 
towards  that  for  conventional  threshold  levels.  The  psychometric  function 
is  then  simply  the  cumulative  (error)  function  derived  from  the  second 
noise  process  (equation  7  and  Annex).  Therefore 

0/  (q"Ki . K2^"° • 28  (9) 

This  value  of  o/Cq  compares  favourably  with  independent  and  direct 
estimates  produced  by  other  authors.  For  example,  Blackwell  (Ref  42) 
derived  values  in  the  range  0.314  to  0.584  from  data  obtained  by  36 
observers. 

An  example  of  practical  application  of  the  contrast  discrimination  model 
is  now  considered.  As  discussed  in  the  Introduction,  there  is  an  inter¬ 
action  between  spatial  frequency  and  the  number  of  detectable  contrast 
steps  (grey  scale).  This  can  be  easily  appreciated  by  considering  the 
detection  of  a  fine  grating  which  has  a  spatial  frequency  approaching  the 
resolution  limit.  At  such  high  frequencies  only  zero  contrast  and  a 
contrast  of  one  need  to  be  represented.  The  required  number  of  contrast 
steps  is  therefore  some  function  of  the  conventional  threshold  level  Cq. 
However,  contrast  discrimination  AC  is  clearly  not  equal  to  Cq,  irrespec¬ 
tive  of  the  base  contrast  C  as  assumed,  for  example,  by  Roetllng  (Ref  29). 
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PIG  11  Theoretical  calculation*  obtained  using  equation  19  (Annex). 

(a)  Generalised  contrast  sensitivity  functions.  Each  curve  Illus¬ 
trates  the  contrast  difference,  as  a  function  of  spatial  frequency, 
obtained  for  a  fixed  level  of  base  contrast.  Values  of  base  contrast  are 
shown  adjacent  to  the  corresponding  curves.  Data  for  zero  base  contrast, 
le  the  conventional  contrast  threshold  levels,  are  taken  from  Campbell  A 
Robson  (Ref  4). 

(b)  Umber  of  discrlnlnable  contrast  steps  calculated  as  a  function 
of  log. spatial  frequency.  The  nuaber  of  steps  Is  also  represented  on  a 
base  2  logarithmic  scale. 

(c)  As  (b)  except  chat  spatial  frequency  Is  represented  on  a  linear 
scale. 
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On  the  basis  of  the  relationship  between  AC  and  C  given  in  Figure  8c,  the 
contrast  sensitivity  function  for  sinusoidal  gratings  can  be  generalised 
to  give  contrast  difference  functions.  The  contrast  threshold  (ie  the 
conventional  threshold  level  Cq)  is  not  determined  as  a  function  of 
spatial  frequency;  instead  the  contrast  is  fixed  at  some  level  C,  and 
values  of  AC  are  determined  as  a  function  of  spatial  frequency.  Different 
values  of  C  then  produce  different  curves  of  AC  against  frequency. 
Conventional  contrast  threshold  values  were  obtained  by  fitting  a  smooth 
curve  through  data  given  by  Campbell  &  Robson  (Ref  4).  The  curve  was 
used  to  produce  the  generalised  contrast  sensitivity  functions  shown  in 
Figure  11a.  The  cor r res ponding  required  number  of  contrast  steps  at 
different  spatial  frequencies  is  plotted  in  Figure  lib  on  a  log  frequency 
scale  and  in  Figure  11c  on  a  linear  frequency  scale.  Each  value  was 
calculated  by  counting  the  number  of  contrast  discrimination  steps  from 
the  conventional  threshold  level  up  to  a  contrast  of  1.0.  The  values  are 
also  shown  on  a  base  2  logarithmic  scale  in  order  to  illustrate  the 
number  of  hits  required  at  each  spatial  frequency.  Note  in  particular 
that  the  number  of  bit  levels  is  greatly  reduced  at  medium  to  high  spatial 
frequencies.  The  curves  shown  in  Figure  11a  illustrate  that  care  must  be 
exercised  when  applying  contrast  sensitivity  functions  (conventional 
threshold  values,  0»0)  to  suprathreshold  situations.  For  example,  the 
process  of  recognition  of  real  objects  often  requires  the  detection  of 
features  within  the  objects  which  are  themselvea  superimposed  on  other, 
suprathreshold  features.  In  addition  calculations  of  the  number  of 
detectable  contrast  steps  (Figs  llb&c)  suggest  that  novel  and  efficient 
ways  of  image  coding  may  be  possible  which  produce  no  detectable  changes 
in  the  image  when  subsequently  decoded.  Finally,  improved  signal  transfer 
functions  can  be  devised  for  imaging  systems  whlchh  use  some  intermediate 
digital  storage. 

5  OONCUJSIONS 

a)  The  Just-noticeable-difference  (JND)  in  contrast  of 

targets,  localised  in  both  space  and  spatial  frequency,  is 
approximately  8X  when  the  target  itself  is  detectable. 

b)  The  function  describing  the  variation  of  JND  with  contrast 
level  is  approximately  invariant  for  wide  ranges  of  experimental 
conditions,  provided  that  both  variables  are  scaled  with  respect 
to  che  conventional  target  threshold  level. 

c)  The  results  can  be  applied  immediately  to  the  specification 
of  Just  detectable  changes  in  imaging  system  MTF  and  to  the 
required  accuracy  in  the  measurement  of  MTF. 

d)  The  experimental  results  and  theoretical  description 

provide  a  firm  basis  for  the  construction  of  visual  models 
applicable  to  the  prediction  of  tho  detectability  of  changes  in 
imaging  system  parameters  or  desig.is,  and  to  the  detection  of 
targets  or  features  within  target?,  especially  when  presented  on 
structured  background  fields. 
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THEORETICAL  HDDEL  OF  CONTRAST  DISCRIMINATION 


ANNEX  A 


Assume  that  the  noise  processes  generste  Gaussian  probability  density 
distributions  when  referred  back,  to  Input  contrast  levels.  Thus,  an 
input  of  contrast  Cq  produces  an  output  which  at  some  time  and  with 
probability  P(C';Cn,c)dC'  Is  identical  to  an  Input  of  contrast  C'  where: 

P(C* 5(^,0)  •  (l/o/5w).expt-(C,-Cn)2/202]  (10) 

Two  targets,  contrasts  C  and  C+AC,  generate  the  two  probability  density 
distributions  shown  In  Figure  10b.  The  dashed  line  at  contrast  Cq 
represents  the  threshold  device.  Responses  which  correspond  to  input 
contrasts  below  Cq  are  not  transmitted  through  the  threshold  device.  The 
forced-choice  experiments  determine  the  probability  that  a  target  of 
contrast,  OAC,  Is  correctly  assigned  as  being  of  higher  contrast  than 
one  of  contrast  C  for  AC  >  0.  This  probability  function  Is  composed  of 
three  parts.  Firstly,  the  following  equation  gives  the  probability 
that  both  targets  are  detected  and  also  that  a  correct  choice  of  target 
is  made. 


?!  -  rP(C' ;C,oc)  /*P (  C" ; C+AC , 0q+^c ) • dC" . dC '  (11) 


Secondly,  the  probability  P2  that  the  target  of  contrast  C  is  not  detected 
but  that  of  contrast,  OAC,  Is  detected  la: 


P 2  m  [  °P(C'  ;C,oc).dC'  .£*P(C'  jC+AC.o^^c) .dC'  (12) 


This  again  corresponds  to  the  probability  of  a  correct  choice.  Thirdly, 
the  probability  P3  that  both  targets  are  not  detected  and  a  correct  choice 
is  made  by  guessing  la: 


V  V* 

P3  -  if  °P(C' ;C,Og)adC. /  0P(C’ ;C+AC,ac+AC).dC' 


(13) 


For  all  conditions  the  probability  P4*  of  a  correct  choice  is  the  sum 


P+  -  P  i+P  2+p  3 


Using  erf(x)  •  (2//r)./exp(-y2).dy 

0 


(14) 

(15) 


"  n-Cn 

than  /  P(C*  jC^oKdC'  -  *.ll-erf( - >’ 

x  /2a 


(16) 
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and  equation  14  can  be  rearranged  to  yield: 

F4,  -  ’+l.ll4«rf(50CC)l  .  1 1-arf ( CnrO-ACM 

^°OfAC 


„  X  •  — - —  •  /  «xpl“(C'-C)2/2oc2)  .  erf  / C*-Q-AC^dC' 
2  /2ioq  Go  Of  AC 

Scale  all  contrast  teres  with  respect  to  by  setting: 

T  -  C/Cq 
AT  -  A<V  Cq 
ST  -  op/ Co 

and  S^x  -  a^Q/ Cq 
Then: 

p+  -  i  +  i  .  [lferf/l-T  xl  .  f  1-erf  /  1-T-/.T  vl 

2  »  <T1S?  (77w 

m 

-  JL  .  — i —  .  /  exp[-(x-T)2/2ST2]  .  erf/X-T-AT  \.dx 
2  /5xSt  1  ^sTfAT 

Prom  equation  8: 

ST  -  Kj.[K2+T2]^ 


(17) 


(18) 


(19) 


and 


(20) 


STfAT  *  ^.{^(TfAT)2!* 

The  probability  of  correct  discrimination  P*  Is  dependent  only  upon  values 
of  contrast  which  are  scaled  with  respect  to  Cq,  the  conventional  threshold 
value.  If  K}  and  K2  are  Independent  of  Cq,  then  the  curve  of  log  contrast 
difference  log. AC  versus  log  base  contrast  log.C  for  a  given  probability 
of  discrimination  has  a  fixed  shape  and  can  be  adjusted  for  targets  having 
different  values  of  Cq  by  sliding  along  a  45deg  line  as  described  In 
Section  3. 

The  theoretical  curve  of  AC  against  C  shown  in  Figure  8c ,  was  produced  by 
solving  equation  19  numerically  for  a  given  value  of  T  and  adjusting  AT 
In  a  systematic  manner  until  F4-  equalled  0.707.  Theoretical  probability 
of  discrimination  values  shown  in  Figure  9  were  generated  for  a  given 
value  of  T  by  a  straightforward  numerical  solution  of  equation  19  with 
different  values  of  AT,  each  value  yielding  a  corresponding  probability 
F*. 
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Appropriate  values  for  the  constants  Kj  and  K2  In  equation  20  were  selected 
separately.  For  large  T,  S?  is  approximately  equal  to  KjT.  Therefore, 
the  constant  K}  can  be  obtained  from  equation  19  with  the  condition  that 
AC/C,  and  hence  AT/T,  la  equal  to  0.08,  the  mean  value  of  the  contrast 
Webor-Fechner  fraction  (Fig  7).  The  value  of  K2  was  then  Increased  from 
sero  for  C/Cq-T*1  ,  as  shown  by  the  asterisks  In  Figure  8c,  until  satisfac¬ 
tory  agreement  was  obtained  with  the  experimental  data. 
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